Comment on "Diffusion of Ionic Particles in Charged Disordered Media" 
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In a recent Letter, Mehrabi and Sahimi discuss motion 
of ions in charged disordered media, presenting a vari- 
ety of results obtained by Monte Carlo simulation on a 
lattice pi. Their observations of mean square displace- 
ments iP(<) suggest that this model exhibits anomalous 
diffusion in three dimensions: 
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They observe the same behavior in one and two dimen- 
sions, but do not present results for 5. Mehrabi and 
Sahimi also make the physically-surprising claim that a 
suitably-defined "short-time diffusivity" can actually in- 
crease with increasing disorder strength (see Fig. 3 of 
Ref. |l|). Exact bounds, renormalization group calcula- 
tions, and previous numerical simulations are inconsis- 
tent with these results in three dimensions. 

At low concentrations of mobile ions, the Green func- 
tion for a diffusing ion should obey the diffusion equation: 


dcv{r,t) 2 

— = DoV Cy 

at 


l3DoV-[cyVv{r)] 


(2) 


Here c„ is the Green function of a single ion in a given 
realization of the quenched random potential v, Dq is 
the "bare," short-time diffusivity, and (3 is the inverse 
temperature. The mean square displacement is given by 
Ry{t) = J d'^r\r\'^Cv{r,t). The observable mean square 
displacement is given by an average over all realizations 
of the disorder: R^{t) = (i?J(t)). The effective diffusion 
diffusion coefficient is defined in d dimensions by D = 

Mehrabi and Sahimi model the disorder by a quenched 
Gaussian random potential field. The statistics of this 
potential field are chosen so that they obey bulk charge 
neutrality: x„^(k) = 7/[fc^(fc^ -|-k^)]. Here the potential- 
potential correlation function is x„^(r) ~ {v{0)v{r)), k is 
an inverse correlation length, and 7 is a measure of the 
density of defects. The Fourier transform in d dimensions 
is given by /(k) ~ J d'^rf{r) exp(zk • r). 

The single-ion, random diffusion model is a well- 
studied one in statistical physics, and a variety of exact 
results are known. First, there is an exact bound for the 
diffusivity in this system in any dimension B : 
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Calculating this bound in three dimensions, one finds 
D/Dq > exp[— /3^7/(47rK)]. This result implies that the 
motion is diffusive in three dimensions, i.e. D > 0. The 
motion is also diffusive at finite ion concentrations, since 
the dynamical exponent is 2 [p[. Therefore, the motion 


should be asymptotically diffusive in three dimensions. 
Indeed, previous careful simulations by Dean, Drum- 
mond, and Horgan on related models have confirmed the 
bound [0. Moreover, these simulations have shown that 
Deem and Chandler's single-ion prediction [^ 
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is accurate to at least moderate disorder strengths. In 
fact, this equation is correct to second order in /3^Xt,t,(0) 
in all dimensions and is exact in one dimension. Note 
that, as expected physically, the diffusion constant de- 
creases with increasing disorder strength. 

The situation is more interesting in two dimensions, 
where anomalous diffusion can occur [the bound in Eq. 
(y) vanishes]. Indeed, field-theoretic treatments have 
shown that the exponent in Eq. (|l|) is continuously vari- 
able and is given exactly by i5 = 1/[1 -I- 87rK^/(/3^7)] [||. 
This scaling has been confirmed by numerical simulations 
W. At finite ion concentrations, the anomalous diffusion 
persists at high temperature 0] , although the mobile ions 
may partially screen the disorder. A Kosterlitz-Thouless 
transition can occur at low temperature P]. 
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